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ABSTRACT 


Selective problem areas relating to computer based ESM 
system software development are discussed. Digital simula- 
tion is presented as a practical approach to solving software 
development problems. An algorithm to simulate a hard-wired 
EW signal preprocessor is developed. An application oriented 
high level language version of the algorithm is presented. 

A FORTRAN program which generates realistic data for 
Signal analysis is described. Using the simulated data, an 
instruction level simulation program, and a language compiler, 


the algorithm is tested using an IBM 360 computer. 
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I. INTRODUCTION 


A. STATEMENT OF ‘THE PROBLEM 

Electronic Warfare Support Measures (ESM) is that 
Givision of Electronic Warfare (EW) concerned with exploita- 
tion of enemy use of the electromagnetic spectrum. ESM 
includes actions taken to search for, intercept, locate, and 
analyze electromagnetic radiations to support tactical 
military operations. 

The continuously increasing density and complexity of 
today's electromagnetic environment, created by friendly as 
well as hostile forces, has fostered the requirement for more 
effective performance and greater sophistication in operational 
ESM systems. Current receiver technology has kept pace with 
this escalation, and signal intercept is virtually assured. 
However, the resulting increased data rates, existing and 
forecast manpower limitations, and limited reaction time 
available to the operational commander have imposed a require- 
ment for faster and more accurate Signal analysis, classifica- 
tion, and interpretation. Currently, major efforts are being 
put forth to achieve an effective ESM system design incorpo- 
rating high speed digital computers for the purposes of 
automatic signal analysis and system control. 

To perform effectively in present and future environments, 
an automated ESM system should exhibit the following 


characteristics: 
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. ‘Broad frequency range. 


. Simultaneous coverage of entire frequency range. 


1 
2 
a eah probability of intercept. 
4. Low false alarm rate, 

5. Capability to recognize pulse or radar-like signals. 
6 


. Capability to determine parameters associated with 
radar pulse trains. 


fe. Capability. to classify szgnals: according to type and 
source. 


See oireceion Linding ‘capability . 
9. Capability to handle data rates. 
ZO, | Capability to detect signals at ranges to radar horizon. 


The successful design and development of a system exhibiting 
these characteristics rest upon three equally important 
elements: system hardware, system software, and the man- 


machine interface. 
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l. ESM Hardware/Svstem Description 


A discussion of ESM hardware as a specific problem 
area is beyond the scope of the thesis. A system concept 
compatible with the state of the art approach is presented 
here as a basic foundation for this study and future efforts 
in the ESM arena. 

Conceptually, an airborne ESM system with the charac- 
teristics previously described could be implemented utilizing 
a combination of Instantaneous Frequency Measurement (IFM) 
and Superhet receivers. The IFM’receivers, operating in 
conjunction with a high-gain pencil beam antenna, would cover 
the frequency bands of interest and function as the primary 
acquisition receivers of the system. In addition to contrib- 


uting to increased system sensitivity, the high-gain pencil 
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beam antenna functions to reduce the instantaneous data rate 
and, in conjunction with azimuth centroiding, provides an 
accurate means of determining Signal direction of arrival. 

Additional data rate reduction is accomplished by 
routing each IFM output through a Pulse Train Separator (PTS), 
which functions as a signal preprocessor to separate stable 
pulse trains from the received data. Parameter measurements 
are performed on those stable trains separated and a computer 
word is generated to tag the associated emitter. The required 
capacity of the PTS in number cof pulse trains separated is 
considered paramount to system operational effectiveness. 

A computer word is generated for each pulse associated 
with non-stable pulse trains and those stable pulse trains 
exceeding the capacity of the PTS. These data are routed 
through a residual channel of the PTS to a high-speed digital 
computer for deinterleavino and parameter measurement. The 
parameters associated with each signal are then passed to a 
main frame computer for classification, threat analysis, and 
tactical environment update. 

The IFM system employs an omni antenna for side lobe 
Suppression to prevent signal poke~-through in the receiving 
antenna minor lobes. The IFM sensitivity and antenna gain 
provide a capability for radar detection out to the radar 
horizon for aircraft altitudes in the neighborhood of 
30,000 feet. 

Priority signals, unclassifiable signals, and signais 


too complex for automatic analysis are automatically assigned 


10 


® a ¢< 
— i 
“e ae. 


= . 
i. a 
twhrenm tastes sf2° esutead® oy : : 
@ ry = 
img ,pnlblesaiss dgumian dete” viva s 
_ a a 


‘sot ity LEAS vatininsotel 3a A 


Ww 


ion -atgorbe «20% ste GQenptsLaag) 
—o . cpeeds stoqdire pOET de 
coger? totete 6s aa0idg 

Poh. wvdeasos =ig mon? e& 

~-4 ofds¢= arodt oo tewtale 

' p< u'4 gaz of hes? -onep 

Hq 3 (odene al 9 ey. Io VW 


21242 82 Je Tee 
a 


teres 0) wa gigdence & 


, 9 


i” eclece ali olde Se—0e 
-P. Sees ors’ 
io le ;eeky BBicows 
‘2ivseliteasaep -eF 
ST (gdCee%e RISE 
o's 30 SA9Neriscg ane ¥ 
drweret. . iow 
ie oppo Pa WTS ae 
'a#eqoso of efika t 
(ir ent redoi, Torntin 


Y ysttirdeqes 8 


900 (gs 47942sale Wi 


“Wad ws 4 


any Ps, 
"nA ee 9 volun 


a 
5 


to digitally controlled superhet receivers for refined 
analysis. The superhet receivers cover the frequency bands 
of interest and operate in conjunction with an omni or low- 
gain broad beam antenna. 

2. Man-Machine Interface 

The human operator element of an automated ESM system, 
though normally an area addressed by human factors engineers, 
impacts heavily on the software requirements and will be 
discussed here in that regard. 

Automation does not eliminate the requirement for 
human operators, but it does significantly affect the role 
the operator plays in the ESM system and it may seriously 
retard the operator's proficiency and ability to function 
efficiently as a back up in the event of a system failure. 
Thus, in deciding which system functions should be automated 
and which should be assigned to the human operator, it is 
imperative that the operator be assigned those functions that 
maximize the use of his capability to exercise judgment, 
improvise, reason inductively, and profit from experience. 

The human operator must be able to control and direct 
the operation of the system. To do so, he must be provided 
relevant information by the system at the proper time, in the 
proper format, and on the proper displays. To obtain efficient 
and reliable software, these requirements must be clearly and 
accurately defined in the planning stage of the ESM system 


development. 
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3. ESM Software 

It is felt that software represents the most trouble- 
some problem in the design and development of an automatic 
ESM system. AsSuming cost and risk vary Girectly with 
complexity, ESM software will account for a substantial 
percentage of total system costs and be the riskiest part of 
the system development. Since system performance hinges on 
software performance, implementation of an effective software 
test plan would serve a threefold purpose of cost reduction, 
risk reduction, and increased system performance. 

The ability to successfully implement such a test 
plan requires a detailed description of the systems functions 
and a software modular building process. Each software 
module or set of modules is then integrated into the system 
and tested against sets or subsets of system functions 


according to prescribed test plans. 


Be SCOPE OF THE THESIS 

Simulation, via computer software, is one technique that 
is useful in software development and verification cr testing. 
Realistic data, created by successful simulation of the ESM 
system hardware operating in a dynamic environment, utilized 
as inputs to a high-speed digital computer would provide the 
framework for implementation, testing, and optimization of 
ESM signal processing algorithms. The inherent ability to 
control the density of the electromagnetic environment through 


Simulation would provide the capability to analyze worst case 
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situations, thereby enabling system designers to estimate 
computer capacity and reserve capacity requirements as well 
as data rate reduction requirements. 

Given these simulation capabilities, an instruction-level 
simulation of the target computer provides a powerful means 
of debugging software and algorithm testing for accuracy and 
speed of operation. Applied Technology Division of Itek 
Corporation, Sunnyvale, California, has provided the Naval 
Postgraduate School with the programming support system for 
their ATAC computer. Instruction level simulation is only 
a part of this software package which includes the compiler 
programs for Programming Language for ATAC (PL/ATAC), a high 
level language developed for ATAC applications programming. 
Acquisition of the ATAC Programming Support System (APSS) 
gives impetus to the utilization of simulation as a technique 
for ESM software development and testing at the Naval Post- 
graduate School. 

Total hardware simulation of the ESM system previously 
described is beyond the scope of the thesis. Total system 
simulation does, however, lend itself to a software "build" 
process and as a first effort of the process, three specific 
objectives of this work are as follows: 

1. Provide a digital simulation of an IFM receiver operating 


in conjunction with a high-gain, pencil beam DF antenna 
in a dynamic radar environment. 


2. Provide a digital simulation of a hard-wiredpulse train 
deinterleaver fed by an IFM receiver. 


3. Describe and demonstrate the capabilities of APSS, 
particularly, the PL/ATAC programming language. 
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The overall objective of this work is to provide a building 
block as a basis for follow on study of ESM system scftware 


at the Naval Postgraduate School. 
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II. COMPUTER SYSTEM FOR ESM SOFTWARE DEVELOPMENT 


—— - 


A. GENERAL DESCRIPTION 

APSS was written for Applied Technology's Airborne 
Computer (ATAC). ATAC, Specifically designed for Electronic 
Warfare applications, is well suited for ESM software design, 
development and testing. A number cof features of this 16-bit, 
fixed point, two's complement computer permit the establishment 
of a meaningful baseline for determining the relative merit of 
algorithms applicable to computer based ESM systems. Special 
features of interest are: 


1. 437 nanosecond major cycle time. 


2. Sixteen general registers which may be used as 
accumulators or as index registers. 


3. Microprogremnable processor, permitting optimization of 
the instruction set for specific applications. 


4, Memory types can be mixed in any combination of low power 
Read Only Memory (ROM), fast bipolar Random Access 
Memory (RAM), and slower core memory devices. 


5. Memory expandable to 65K words. 


6. Four basic instruction addressing modes: 
register, immediate, direct, and direct indexed. 


7. Three suppiemental addressing modes: 
register indexed, immediate short, and direct register. 


8. Provisions for up to 16 Direct Memory Access (DMA) 
channels. 


9. Provisions for up to four simultaneously active DMA 
channels. 


10. Twelve hardware priority interrupts. 


11. Modular, compact construction. 
Additional detailed information, including the instruction 
set, may be obtained from ATAC, Applied Technology Airborne 


Computer, Vow. L,, Principles of Operation. 
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B. ATAC PROGRAMMING SUPPORT SYSTEM (APSS) 

This software support system allows for the preparation 
and debugging of programs via an IBM 360/370 computer system. 
APSS, through simulation of the ATAC computer and associated 
peripheral devices, permits the applications programmer to 
optimize programs for execution speed and memory utilization. 
The APSS system is written in a combination of PL/1, FORTRAN 
IV, and IBM assembler language and includes the programs 
listed below. 


APSS Monitor 

ATAC Instruction Assembler 

ATAC Instruction-Level Simulator 
Relocating and Linking Loader 
Microcode Assembler 
Microcode-Level Simulator 
PL/ATAC Compiler 

Library Control System 


Subroutine Library 
Arrangements were made with Applied Technology to use portions 
SGeAPsSs and, consequently, the ATAC. Instruction-Level Simulation 
Program and the PL/ATAC Compiler are currently operational on 
the IBM 360/67 at the Naval Postgraduate School. Information 
concerning the implementation of these programs and instruc- 
tions pertaining to the use of the programs are contained 
in Appendix A. Additional APSS information is available in 
ATAC, Applied Technology Airborne Computer, Volume 2, ATAC 
Programming Support System (APSS). The documents in this 
volume provide a detailed description of how to set up and use 


the APSS system. Those documents applicable to this work are: 
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1. Simulation of instruction execution. 


2. Tracing option of the 16 general registers. The contents 
of each register can be printed after the execution of 
an instruction. 


3. Accumulated program execution time printed with each 
line of trace. 


4, Instruction execution histograms by frequency of 
execution or time of execution can be created and printed. 


5. Partiai and total simulated memory dump options. 

6. Initialization of program data arrays allowing data to 
be read into virtual program memory prior to program 
Simulation. 

The APSS Instruction-Level simulator program includes a large 
repertoire of error messages to assist the applications 


programmer in debugging job management and ATAC instruction 


uSage errors. 


C. PL/ATAC 

PL/ATAC, Programming Language for ATAC, is.a high level, 
block structured language designed specifically for the ATAC 
computer. The programming structure and numerous capabilities 
cf the lanquage are demonstrated in Appendix C via the Pulse 
Train Separator simulation program written in PL/ATAC. In 
addition to hardware simulation, a specific objective in 
writing this program was to provide programming examples of 
the PL/ATAC language to assist the interested programmer in 
learning PL/ATAC. The PL/ATAC manual (Volume 2, Section II, 
ATAC) is an excellent reference but may present some difficulty 
if uSed as a learning text. 

PL/ATAC is a powerful language providing the programmer 


with programming control approaching that of assembly language 
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without compromising the high level of the compiler language. 
Additionally, ATAC assembly language may be intermixed with 
PL/ATAC increaSing programming flexibility and control with 
no serious degradation of program readability. 

An analysis to determine the effectiveness of applications 
programming in the PL/ATAC language versus ATAC assembly 
language was conducted by the author utilizing an assembler 
optimized ATAC benchmark program. The results of this analysis 
indicated a memory utilization increase by a factor of 
approximately 2.3 times the ATAC assembly language requirement 
and an execution time increase by approximately 2.4 times the 
ATAC assembiy language execution time. 

The PL/ATAC programmer loses some flexibility, as compared 
to the assembly language programmer, due to the limitation in 
the number of registers available (a maximum of five) to him. 
Also, the PL/ATAC compiler currently installed at the Naval 
Postgraduate School does not have a double precision capability. 
Double precision calculations, if required, must be accomplished 
in assembly language resulting in increased core space and 


execution time required. 
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IIIT. ESM ENVIRONMENT/IFM SIMULATION 


Investigation revealed the existence of a computer program 
[Ref. 1] designed to simulate an IFM receiver operating in 
conjunction with an omni antenna in a dynamic redar environ- 
ment. This program was modified to incorporate a rotating 
DF antenna system. Routines were added to generate unstable 
pulse trains in the form of pulse 'jitter' and/or 'stagger'. 
Documentation was added to facilitate program readability. 

The DATA Simulation Program generates simulated pulsed 
radar Signal level data arriving at an IFM receiver via a 
rotating antenna system. Current capacity aliows for simula- 
£10n of up te ten emitter platforms with up to five emitters 
per platform. This capacity may be increased by increasing 
the array dimensions of the parameters associated with each 
emitter-platform combination. Emitter platform and IFM 
Platform initial position information in X-Y grid coordinates, 
course in degrees clockwise from North, and speed in knots 
are set by the program input data. 

“Emitter characteristics, set by the program input data, 
consist of transmitter frequency in MEz, pulse repetition 
interval (PRI) in milliseconds, pulse width (PW) in micro- 
seconds, transmitter power in kilowatts, antenna gain in db, 
antenna beamwidth in degrees, antenna scan rate in seconds 
per revolution, and pulse jitter and/or stagger characteristics. 


A sin(x)/x antenna radiation pattern is assumed for ail emitters. 
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IFM receiver characteristics, set by the program input 
data, consist of receiver threshold sensitivity in dbm, upper 
and lower frequency limits of the passband in MHz, minimum 
time measurement capability in microseconds, antenna gain in 
db, antenna scan rate in Seconds per revolution, and antenna 
beamwidth in degrees. 

Additional program inputs include simulation start and 
stop times in seconds, frequency bin size of the IFM digital 
processor, and random number generator calling arguments. 
Complete input data description is depicted in the program 
listing following the description of the program variables. 

A complete program listing is contained in Appendix B. 

Utilizing the input data described, the Data Simulation 
Program computes the position, course, speed, and range of 
each emitter platform relative to the IFM platform. The 
closest point of approach (CPA) and distance to the CPA is 
computed for each emitter platform to determine if and when 
a platform will be within the range of detection of the IFM. 
Those emitters out of the detection range during the simulation 
are listed as such. The IFM antenna's and all emitter 
antennas' intial pointing directions are set by a random 
process utilizing a uniform random number generator. The time 
at which each emitter platform enters the capture area of the 
IFM antenna's main beam is calculated and used as a basis for 
computing a random transmit time for the emitters associated 


with that platform. 
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Emitter maximum effective radiated power (P._) is then 


E 
computed utilizing the mathematical model given by 


és 2 2 
Pa = P. G, A“ /16T 


where 
Py = Transmitted Power, 
G. = Emitter Antenna Gain, 


and > = Wavelength of Signals. 


Emitter maximum range of detection (R_) is computed by: 


1 
es bs 
Ra (PE, G/T) 


where 
G.. = IFM Antenna Gain, 
and Th = IFM Receiver Threshold Sensitivity. 


For each emitter within detection range, the minimum antenna 
gain (San? required to break the IFM threshold sensitivity 


is then computed by 


e Zz 
Simin i Ty, /G, oe 
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Using the value of emitter transmit time previously calculated, 
the emitter antenna pointing direction is calculated. Emitter 
antenna gain as a function of the relative angle between the 
emitter antenna and the IFM antenna is then calculated and 
compared with the minimum antenna gain required to break the 
IFM threshold sensitivity. If the signal level is sufficient 
to break the IFM threshold, the time of arrival (TOA) of the 
pulse at the IFM receiver is calculated. 

If the pulse does not break the IFM threshold, the emitter 
transmit time is incremented by the emitter PRI, the antenna 
position is updated, and the antenna gain calculation process 
is repeated until a Signal level of sufficient strength is 
transmitted. Should no pulse be received in the scan time 
of one beamwidth of the IFM antenna, the subsequent scan of 
the receiver antenna is checked and the process continues 
until a pulse is received or until simulation stop time. 

The time of arrival of the first pulse from each emitter 
is stored in an array and, upon completion of this process 
for all emitters, the array is searched for the lowest time 
of arrival. The receiver antenna bearing at this TOA is 
calculated and provided as output data along with the parameters 
of the associated emitter. The emitter platform position is 
updated and the time of arrival of the next pulse from this 
emitter @s calculated and stored in the array. The array is 
then searched for the next lowest time of arrival and the 


process repeated until simulation stop time. 
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The output of the Data Simulation Program consists of 
parameters associated with each pulse received by the IFM 
receiver in the proper time sequence. These parameters 
include the frequency, pulse width, bearing and time of 
acriVval. 

This data, as prepared for the ATAC simulation programs, 
consists of five integers for each pulse. These integers 
represent the following: 

1. Emitter bearing. An integer in the range of 0-359. 
This word represents the IFM antenna bearing at the 
time of arrival of the pulse. 

2. Frequency bin. The bins are numbered in ascending 
order of frequency from 1-128. This range may be 
altered as desired by the programmer. 

3. Pulse width. The pulse width word is an integer 
in the range 0-65,535. The time scaling factor is 
arbitrary. Current convention is that the least 
significant bit of the pulse width word corresponds 
to .1 microseconds. This allows for representation 
of pulse widths from 0-6.5 milliseconds to within 
0.05 microseconds. 

4. Time of arrival. Time of arrival is expressed as 
a double precision (32 bit) integer. Word 4 is 
the least significant half of this integer and holds 
time values up to 0.84 seconds with the least signi- 


ficant bit corresponding to 25.6 microseconds. 
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ee Dame. Orlarrival, ~The upper half of the time word 
holds time values up to 15.3 hours. 
The Data Simulation Program produces these outputs in 
printed format and in a punched card format suitable for 
acceptance by the ATAC simulation program. A sample of 
the program output is contained in Appendix B following the 
program listing. The program flow chart is illustrated in 


Figures 1 and la. 
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INITIALIZE CONSTANTS AND SCALE FACTORS 
READ CALLING ARGUMENTS FOR RANDOM NUMBER GENERATORS, 


SIMULATION START AND STOP TIME, IFM CHARACTERISTICS, 
AND IFM PLATFORM COURSE, SPEED, AND POSITION. 


INITIALIZE IPM ANTENNA POINTING DIRECTION 


| OUTPUT I IFM CHARACTERISTICS | 


7 READ EMITTER PLATFORM INFORMATION 


<CPLATFORM ALTITUDE LT ZERO +) 
i] wana AEE TY A Re 


' No 
uy 
OUTPUT EMIT EMITTER PLATFORM INFORM RMATION 


COMPUTE EMITTER PLATFORM COURSE, SPEED, 
AND RANGE RELATIVE TO IFM PLATFORM z 


ease! 2 Ser mera ose RE 


COMPUTE BEARING, CPA, AND DISTANCE TO 
CPA OF EMITTER PLATFORM 


ee ah hal 
COMPUTE TIME WINDOW IN WHICH THIS PLATFORM 
WILL BE IN RECEIVER ANTENNA BEAMWIDTH 


READ EMITTER CHARACTERISTICS, CONVERT TO INTERNAL 
FORM, AND COMPUTE MAXIMUM EFFECTIVE RADIATED POWER 


| ~ OUTPUT EMITTER CHARACTERISTICS — 
eas ae 


COMPUTE MAXIMUM RANGE AT WHICH IFM CAN DETECT EMITTER | 


| COMPUTE TIME OF ARRIVAL OF FIRST PULSE FROM EMITTER — 


aN salar cid 
IS THIS THE LAST EMITTER ON THIS | 


Yes 


FIGURE 1. Data Simulation Program (Simplified Flowchart) 
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| FIND LOWEST TOA AND OUTPUT PULSE DATA FOR THIS EMITTER | 


COMPUTE TOA OF NEXT PULSE FROM THIS EMITTER 
AND UPDATE EMITTER PLATFORM POSITION 


No a abla x 
~ SIMULATION STOP TIME 


—e 


Yes 


STOP 


FIGURE la. Data Simulation Program (Simplified Flowchart) 
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IV. THE PULSE TRAIN SEPARATOR 


A. CONCEPT OF OPERATION 

Data received by the IFM are routed to the Pulse Train 
Separator (PTS). If a pulse train has a stable pulse repe- 
Euerom £Eequency (PRE) the PTS locks on the train and tracks 
the train for a period of time determined by the receiver 
antenna scan rate and beamwidth. The signal's pulse width 
(PS) and pulSe repetition interval (PRI) are measured and 
Gigitized along with azimuth, time first seen and frequency. 
These five parameters are used to tag the associated emitter 
and are routed to the main frame computer for emitter classi- 
fication. Direction of arrival is computed by azimuth cen- 
troiding based on the receiver antenna bearing at the time 
Of arrival of the first and last pulse of the signal. The 
possibility exists that, due to sin(x)/x characteristics 
of emitter antennas, centroiding will be performed on ragged 
Signal inputs. This is not considered degrading to system 
performance as meaningful direction of arrival information 
wil, still: result. 

An inherent characteristic of the Pulse Train Separator 
is that those stable trains with the lower values of PRI will 
be separated first. This feature of the PTS is also desirable 
as the data rate reduction role of the device is enhanced. 
The required capacity of the device in this role of data rate 


reduction has not been eStablished and is considered a most 
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important factor in the overall ESM system design. Failure 
to provide the proper level of preprocessing could result 

in system software and hardware overloading and, consequently, 
degraded system performance. 

Those stable signals exceeding the capacity of the Pulse 
Train Separator and non-stable signals such as frequency 
agile signals, jittered pulse trains, staggered pulse trains, 
and other exotic signals are routed through a residual signal 
channel of the pulse train separator to the high speed mini 
computer for parameter meaSurement. For each of these pulses 
received, a computer word will be generated for frequency, 
time of arrival, pulse width, and direction of arrival. Pulse 
deinterleaving and such parameter meaSurements as pulse 
repetition interval and pulse jitter and stagger character- 
istics will be performed by the mini computer for emitter 


classification. 


B. PTS SIMULATION PROGRAM 

The Pulse Train Separator Simulation Program is designed 
to simulate a hardware preprocessing system with a capability 
of separating wp to a maximum of ten stable pulse trains from 
the output of an IFM receiver. 

Input data to the program consist of a series of five 
integer words for each pulse generated by the Data Simulation 
Program; direction of arrival in degrees clockwise from North, 
frequency by bin number, double precision time of arrival 


(two words) in seconds. withthe least significant bit 
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corresponding to 25.6 microseconds and pulse width in seconds 
with the least significant bit corresponding to .1 microseconds. 
These data are loaded into a large input buffer prior to the 
start of the program simulation. 

Each pulse within the buffer is compared in frequency 
with previously received pulses within the same receiver 
antenna beamwidth. If no match in frequency occurs the pulse 
is considered a new Signal and a pointer address is stored 
in a New Signal Pointer File. If a frequency match occurs, 
the pulse width of the signals are compared. Upon obtaining 
a match in pulse width the matched pulse's pointer address 
is stored in a five element array within the New Signal 
Pointer File. Upon receipt of five pulses of a pulse train, 
the program calls a local re-entrant procedure (PRICOMP) 
wherein the stability of the pulse train is determined. For 
purposes of this simulation a pulse train is considered 
unstable if the variation in time interval between pulses 
is greater than 51.2 microseconds and/or three or more 
missing pulses are encountered in ragged signals. Upon 
een ana a pulse train to be stable, the pulse repetition 
interval of the signal is calculated and the signal is 
flagged as a stable train. The number of available channels 
in the Pulse Train Separator is then reduced by one and the 
signal is zeroed in the input buffer as each additional pulse 
is received. Parameters associated with this pulse train 
are stored and the signal bearing is updated as each pulse 


is received through incrementation of the input buffer. 
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Upon receipt of the last pulse expected from the pulse 
train, the centroid direction of arrival is calculated and 
the signal information is zeroed in the New Signal Pointer 
File as well as in the first five input buffer lcoations. 
The time of arrival of the last pulse expected is computed 
utilizing the scan time for one receiver antenna bandwidth. 

Non-stable signals are flagged as such and remain in the 
input buffer. Upon completion of incrementation through the 
input buffer, all stable pulses, except those from trains 
that remain due to PTS capacity having been exceeded, have 
been zeroed from the input buffer. Parameters associated 
with these trains are stored in sequence according tc pulse 
repetition interval. 

The pulse information remaining in the input buffer, 
associated with the unstable pulse trains and with pulse 
trains exceeding the PTS capacity, is then compressed within 
the input buffer, thus allowing for simulation of the mini 
computer signal processing algorithm. In the initial design 
concept of this simulation program, these remaining data, 
representing the Pulse Train Separator residual channel, 
would be processed by an ATAC assembly language program 
incorporated as a local procedure within the PTS program 
uSing the code selection head as described in the PL/ATAC 
manual. 

A complete PL/ATAC program listing of the Pulse Train 


Separator is contained in Appendix C. An assembly language 
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program generated by the PL/ATAC compiler is contained in 
Appendix D. A program flowchart is illustrated in Figures 


2 and 2a. 
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FIGURE 2. The Pulse Train Separator Program Flowchart 
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FIGURE 2a. The Pulse Train Separator Program Flowchart 
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V. DISCUSSION OF SIMULATION RESULTS 


A. SUMMARY OF RESULTS 

The Pulse Train Separator as described in Chapter IV 
was Simulated utilizing the ATAC Programming Support System 
and the PL/ATAC Programming Language. The particular hardware 
system simulated was arbitrarily given the capability of 
Geinterleaving and preprocessing up to nine radar pulse trains. 

Input data, provided by the Data Simulation Program, 
consisted of a mixture of stable and non-stable pulse trains 
from eighteen emitters located on six platforms geographically 
positioned so as to be intercepted in one receiver antenna 
beamwidth. The receiver antenna beamwidth was arbitrarily 
set at five degrees, 

Using a receiver antenna scan rate of two Seconds per 
revolution, the Data Simulation Program generated a total of 
365 pulses as IFM output during one antenna main beam sweep 
of the eighteen emitters. These pulses represent approxi-~ 
mately 28 milliseconds of real time data. The Pulse Train 
Separator Simulation Program, in processing nine of the pulse 
trains, removed 277 pulses from the input data and routed 
88 pulses through the residual channel. These 88 pulses 
represent the unstable pulse trains and the lower pulse 
repetition frequency (PRF) Signals. Accurate pulse repetition 
interval calculations as well as accurate azimuth centroiding 


was performed on all preprocessed pulse trains. 
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B. CONCLUSIONS 

The Data Simulation Program provides realistic data to 
the APSS programs. The program produces reascnably complex 
data in the form of jittered pulse trains and staggered 
pulse trains. Pulse trains transmitted at sufficiently low 
power levels are produced with missing pulses due to sin(x) /x 
effects. 

The Pulse Train Separator Simulation Program demonstrated 
the use of the PL/ATAC language and functions as designed to 
preprocess ESM signals. This program is somewhat inefficient 
in its present form for two reasons. First, the limitation 
of five registers allocated to the programmer by the PL/ATAC 
compiler leads to a significant increase in program execution 
time as compared to a program written in ATAC assembly language. 
As the programmer gains experience with the PL/ATAC language, 
this limitation of the compiler may be overcome somewhat by 
programmer ‘tricks’ to gain use of registers. 

The lack of a double precision capability in the current 
version of the PL/ATAC compiler installed on the IBM 360/67 
at the Naval Postgraduate School created a severe handicap 
in writing the Pulse Train Separator Simulation Program and 
also contributes to program inefficiency. Each double 
precision calculation performed requires a program segment cof 
assembly language and the housekeeping involved results in 
increased program execution time. 

Even with these limitations it is felt that the simulation 


programs presented in this work provide the framework for ESM 
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software development and applications programming at the 


Naval Postgraduate School. 


C. RECOMMENDATIONS 
Based upon the experience and results obtained from this 
Study, recommendations for follow-on study in the area of 
ESM software development are as follows: 
1. Continue study of the Pulse Train Separator Program 


and optimize the basic algorithm. 


2. Implement the resulting improved algorithm in ATAC 
assembly language. 


3. Use the Data Simulation Program and optimized hardware 
simulation with APSS to develop and optimize ESM signal 
processing algorithms. 


4. Continue use of the PL/ATAC compiler to develop and 
test hardware simulation and other programs prior to 
implementation in ATAC assembly language. 
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APPENDIX A 
INSTRUCTIONS FOR INSTALLATION AND USE OF THE ATAC 


PROGRAMMING SUPPORT SYSTEM (APSS) ON THE IBM 360/67 
OPERATING SYSTEM AT THE NAVAL POSTGRADUATE SCHOOL 


A. INSTALLATION 
The ATAC Programming Support System is currently loaded 
on the 2321 Data Cell. Data set names and their liocations 


are as follows: 


2. $1615.RDS.LOADLIB = CEL 007 
Zio OL. RDS APSS =e tole, (0:0'5 
3. S615 RDS .PLATAC So Clank MOO 


In the event of failure of the Data Cell these data sets may 
be reloaded from the APSS source tape stored at the W.R. 
Church Computer Center under the name 'ATI.006'. 

File #1 of this tape contains a deScription of the tape. 
File #2 contains examples of the JCL required to create anda 
run APSS. The instructions contained in this file are super- 
seded by the instructions contained here. The following JCL 
is required to list these files: 


//(STANDARD JOB CARD) 

//EXEC PGM-IEBGENER 

//SYSPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT2 DD SYSOUT=A, DCB=(LRECL=80,RECFM=F ,BLKSIZE=80) 
LABEL=(1,SL,,IN),DCB= (BLKSIZE=800,LRECL=80, 
RECFM=FB) , VOL=SER=ATI006 ,DSN=FILEL 

// EXEC PGM=IEBGENER 

// SYSPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT2 DD SYSOUT=A ,DCB= (LRECL=89 , RECFM=F , BLKSIZE=80) 

//SYSUT1L DD UNIT=(2400,,DEFER) ,DISP=(OLD,PASS) , 
LABEL=(1,SL, , IN) ,DCB=(BLKSIZE=800 ,LRECL=80) , 
RECFM=FB , VOL=SER=AT1006 ,DSN=FILE2 
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File #3 of the source tape contains the object module 
library of APSS and is required to build the APSS Instruction 
Level Simulator. File #4 contains the PL/ATAC Compiler object 
module. The following steps are required to build APSS from 
fates’ 3 and 4. 

1. Read files 3 and 4 from tape to disc storage. Special 
arrangements should be made to insure that the disc pack used 
has sufficient unallocated space for this step and the 
following step of this procedure. SPOOL3 has been made 
available in the past and is the disc used in the JCL of this 
example. The load modules from file 3 are moved into a disc 
file named 'S1614.RDS.LOADLIB'. File 4 is moved into a disc 
file named 'S1615.RDS.PLATAC'. The utility control statements 
may be changed to rename the files as they are being loaded 
onto the disc. The JCL statements required are as follows: 


//(STANDARD JOB CARD) 

//STEP1 EXEC PGM=IEHMOVE, REGION=100K 

//SYSPRINT DD SYSOUT=A 

//SYSUT1 DD UNIT=SYSDA,SPACE=(TRK, (40) , ,CONTIG) 

//DA1L DD UNIT=2314 ,DSN=S1615.RDS.LOADLIB, 

J /PSPACHE= (TURK, (5071010) , ;,CONTIG) , 

//DISP=(NEW,KEEP) , VOL=SER=SPOOL3 

//TITAPE DD UNIT=(2400,,DEFER) ,DISP=(NEW,PASS) , 

. LABEL=(3),Sil,,1N) , 

//DCB= (DEN=2,BLKSIZE=800 , LRECL=80 , RECFM=FB) , 
VOL=SER=ATI006 

//SYSIN DD * 

COPY PDS=ATI.APSS.LOADLIB, Los 2314=SPOOL3 ,TODD=DA1, 
FROMDD=TILTAPE, 

FROM=2 400=(ATIO0O6,3), Be Bee eanee HOADIA B 

COPY DSNAME=ATI .APSS.PLATAC , TO=2314=SPOOL3, 
FROMDD=TITAPE, TODD=DA2, * 

FROM=2400=(ATI006,4) , RENAME=S1615.RDS.PLATAC 

//DA2 DD UNIT=2314 ,DSN=S1615.RDS.PLATAC,SPACE=(TRK, 
(20,110) , RLSE) + 

//DISP= (NEW, KEEP) , VOL=SER=SPOOL3 , LABEL=RETPD=360 

f* 
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2. The Library (named 'S1615.RDS.LOADLIB') now on the 


disc must be built into an absolute load module of the entire 


APSS system. The module thus built (named 'S1615.RDS.APSS') 


will be stored permanently on cell 005 of the data cell. This 


file may also be renamed as desired. The following JCL state- 


ments are required. 


//(STANDARD JOB CARD) 

//BUILD EXEC PGM=IEWL, REGION=150K, 

//PARM='OVLY, XREF,LET,LIST,SIZE= (256K,20480) ' 

//SYSPRINT DD SYSOUT=A 

//LIBRARY DD DSN=S1615.RDS.LOADLIB ,UNIT=2314, 
VOL=SER=SPOOL3,DISP=SHR 

//SYSLIB DD DSNAME=3YS1.FORTLIB,DISP=SHR 

//SYSLMOD DD DSNAME=S1615.RDS.APSS, UNIT=2321, 
VOL=SER=CEL005, 

//DISP= (NEW ,KEEP) , LABEL=RETPD=360, 

//SPACE=(TRK, (300,20,2) , RLSE) 

Y7/SGYSUTL DD’ UNIT=SYSDA , SPACE=(TRK, (19,19), ,CONTIG) , 
SEP=SYSLMOD 

//SYSLIM DD* 

INCLUDE LIBRARYEPMUDLN 

CHANGE MSIMEIHESAPDN 

INCLUDE LIBRARY (APSSMON) 

INCLUDE LIBRARY (MSIM4,ASEM5,SIM16A,SIMTR1,SIMI01,GUL) 

INCLUDE LIBRARYEXPLMONN 


OVERLAY Al 

INSERT MSIMUL, *MSIMULA, IHENTRY , IHESAP 

INSERT MINT, IN, OUT 

INSERT IHEDBN, IHEXID 

INSERT IHEBSM, IHECSM 

INSERT IHEBSK, IHEIOX, IHEIOP, IHEDID, IHEDOB 

INSERT IHEBIB, IHEDOA, IHEIOB 

INSERT IHEIOA, IHEOCL 

INSERT IHEBSD, IHEBSF 

INSERT IHEJXS 

INSERT IHEOSD, IHEOST, IHEBSI 

INSERT IHEVPE, IHEDMA, IHEVFB 

INSERT IHEDNC, IHEVFD, IHEVFA, IHEVPD, IHEVPB, IHEVSC 
INSERT IHEVSD, IHEVFE, IHEDCN, IHEUPB 

INSERT IHEVFC, IHEVPE, IHEVPG, IHEVQB, IHEVQC 
INSERTIHEABN, IHEIOD, IHEIOF, IHEPRT, IHEVQA, IHESPRT 
INSERT IHEBEG, IHEERR, IHESIZ. 

INSERT MISET 

OVERLAY Al 

INSERT ASSEM, REWIND, REW72, DSKOUT, CARDIN, DISKIN, 


ERPRT, PRIADD 
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INSERT WRDATA, PRICOM, PRINOR, WRITEX, REFTIT, PREF, 
ERILT 

OVERLAY Al 

INSERT PARMRD, PRESIM 

OVERLAY Al 

INSERT SMLTR, NRMTRM, STRTSM, TPAGE 

INSERT RDCRD, ABNRMT, ARTHFP, TRACE, HGRAM, HGRAMI, 
HGRAMS 

PNSERT IOLNIT ,- ACTIVE, SIMTIM 

ENSERT DEVDTA, ACT, TIME, INT, RAND, DEADT, DEBUG 

OVERLAY A2 

INSERT LEVEL, DMAIOI, DMAIOA, DMAIOD, RIOIO 

iotiho REMACT, DMA, DMATM, REO, RIOTM, RIOINT, INTOLY, 
DMAINT 

INSERT DTRAN, PUTACT, RANDOM 

OVERLAY Al 

INSERT HGPRNT 

OVERLAY Al 

INSERT LAODER 

OVERLAY Al 

INSERT LINK, ENTEXT, SLLH 

OVERLAY Al 

INSERT PLATAC, IOPACK 

OVERLAY SOBJECTEREGIONN 

PolKt OBRGECY, INIT, LIB, RCALPH, RCHEX, RCINT 

INSERT MDATE 

OVERLAY SDUMPEREGIONN 

INSERT SMDUMP, PAGE 

ENTRY MAIN 

NAME APSS 

/* 


3. The PL/ATAC compiler object module on the disc must 


be moved to the data cell. 


//(STANDARD JOB CARD) 

//EXEC PGM=IEBGENER 

//SXYSPRINT DD SYSOUT=A 

//SYSUT1 DD DISP=SHR,UNIT=2314 , VOL=SER=SPOOL3, 
//DSN=S1615.RDS.PLATAC 

//SYSUT2 DD DISP=(NEW, KEEP) , UNIT=2321,VOL=SER=CELO01, 
//DCB= (RECFM=FT , BLKSIZE=7188) , LABEL=RETPD=360, 
//SPACE= (7188, (20,0)) ,DSN=S1615.RDS.PLATAC 

//SXSIN DD * 

GENERATE 


/* 
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4. The Library ('S1615.RDS.LOADLIB') must now be moved 


from the disc to the data cell for permanent storage. This 


library contains modules necessary for Successful execution 


of the PL/ATAC compiler. The required JCL statements to 


accomplish this step are as follows: 


B. 


//(STANDARD JOB CARD) 

//EXEC PGM=IEBCOPY 

//SYSPRINT DD SYSOUT=A 

//SYSUT1 DD DISP=SHR,UNIT=2314,VOL=SER=SPOOL3, 
//DSN=S1615.RDS.LOADLIB 

//SYSUT2 DD DISP=(NEW,KEEP) , UNIT=2321,VOL=SER=CELO07, 
//DSW=S1615 .RDS.LOADLIB,SPACE=(7294,(61,0,14)), 
/ /DCB= (RECFM=UT , BLKSI ZE=7294) , LABEL=RETPD=360 
//SYSUT3 DD UNIT=SYSDA,SPACE=(TRK, (20,5) ) 
V/SXSUT4 DD UNIT=SYSDA,SPACE=(TRK, (20,5)) 
//SYSIN DD * 

COPY OUTDD=SYSUT2,IWDD=SYSUT1 


J* 


USER INSTRUCTION 


The JCL required to compile and execute PL/ATAC programs 


is as follows: 


// (STANDARD JOB CARD) 

//PLATAC EXEC PGM=XPLMON , PARM='ALTER' , REGION=200K 

//STEPLIB DD DSN=S1615.RDS.LOADLIB, UNIT=2321,DISP=SHR, 
VOL=SER=CEL007 

//PROGRAM DD DSN=S1615.RDS.PLATAC, UNIT=2321,VOL=SER=CELOOI, 
DISP=SHR 


-//SYSPRINT DD SYSOUT=A 


//FTO5F001 DD DDNAME=SYSIN 
//FTO6FO01 DD SYSOUT=A 
//FT27F001 DD DSN=&&PLCODE, SPACE= (TRK, (20,10) ,RLSE) , 


Va DCB= (BLKSI ZE=3520 , LRECL=80, RECFM=FB) , DISP= (NEW,PASS, DEL 


Vie UNIT=SYSDA 
//SYSIN DD * 


PL/ATAC Program 
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//APSS EXEC PGM=APSS REGION=200K 
//STEPLIB DD DSN=S1615.RDS.APSS ,UNIT=2321,DISP=(OLD,KEEP) , 
//  ‘OL=SER=CEL005 ,DCB=BLKSIZE=400 
//FT05F001 DD DDNAME=SYSIN 
//FTO6F001 DD SYSOUT=A,DCB= (RECFM=FBA , LRECL=133,BLKSIZE=3325) 
//FTO8FO0L DD UNIT=SYSDA,SPACE= (CYL, 1) ,DCB=BUFNO=1 
//FTO9FOOL DD UNIT=SYSDA , SPACE= (CYL, (7,2)), 
DCB= (RECFM=VBS , BLKSIZE=7180 , LRECL=92 , BUFNO=1) 
//FT1OFO01 DD UNIT=SYSDA,SPACE=(CYL,(7,2)), 
DCB= (RECFM=VBS , BLKSIZE=4204 , LRECL=42 , BUFNO=1) 
//FT20F001 DD UNIT=SYSDA,SPACE=(CYL,(7,2)), 
DCB= (BUFNO=1 , RECFM=VBS., BLKSIZE=2004 , LRECL=500) 
//FT27F001 DD DSN=&&PLCODE, UNIT=SYSDA,DISP=(OLD,DELETE) , 
DCB=BUFNO=1 
//SYSPRINT DD SYSOUT=A , DCB= (RECFM=FBA , LRECL=133 ,BLKSIZE=3325) 
//SYSIN DD * 
SEND 
/* 
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APPENDIX B 


DATA SIMULATION PROGRAM FORTRAN LISTING AND OUTPUT 


A description of program variables is included in the 


Data Simulation Program listing. The input data deck setup 


is also included in the listing. 
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